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ABSTRACT 

NGC 6908, a SO galaxy situated in direction of NGC 6907, was only recently recognized 
as a distinct galaxy, instead of only a part of NGC 6907. We present 21 cm radio synthesis 
observations obtained with the GMRT and optical images and spectroscopy obtained with the 
Gemini North telescope of this pair of interacting galaxies. From the radio observations we 
obtained the velocity field and the H 1 column density map of the whole region containing 
the NGC 6907/8 pair, and by means of the Gemini multi-object spectroscopy we obtained 
high quality photometric images and SA resolution spectra sampling the two galaxies. By 
comparing the rotation curve of NGC 6907 obtained from the two opposite sides around the 
main kinematic axis, we were able to distinguish the normal rotational velocity field from 
the velocity components produced by the interaction between the two galaxies. Taking into 
account the rotational velocity of NGC 6907 and the velocity derived from the absorption lines 
for NGC 6908, we verified that the relative velocity between these systems is lower than 60 
km s~^. The emission lines observed in the direction of NGC 6908, not typical of SO galaxies, 
have the same velocity expected for the NGC 6907 rotation curve. Some emission lines are 
inside their equivalent absorption lines, reinforcing the idea they were not formed in NGC 
6908. Finally, the H 1 profile exhibits details of the interaction, showing three components: 
one for NGC 6908, another for the excited gas in the NGC 6907 disk and a last one for the gas 
with higher relative velocities left behind NGC 6908 by dynamical friction, used to estimate 
the time when the interaction started in (3.4 ± 0.6) x 10^ years ago. 

Key words: galaxies: interactions - galaxies: kinematics and dynamics - galaxies: photom- 
etry - radio lines: galaxies. 



1 INTRODUCTION 

Searching for companions among the nearby barred galaxies from 
the Shapley-Ames Catalog, Garcia-Barreto et al. (2003) identified 
NGC 6908 as a discoidal companion of NGC 6907, 10 kpc away 
from it. They estimated this distance based only on their angular 
separation. But according to these authors, the alignment of the ob- 
jects could be a mere effect of geometrical superposition, since the 
narrow band Ha images of NGC 6907 did not reveal any emission 
in direction of NGC 6908 (Garci'a-Barreto et al. 1996). Surace et 
al. (2004) recognized NGC 6907 and 6908 as a potential interact- 
ing system, but the low resolution of their far-infrared images did 
not allow any further conclusion. 

In a recent paper, Madore et al. (2007) revealed that NGC 



* E-mail: scarano@astro.iag.usp.br 



6908, once classified as a PofG (part of a galaxy), is in fact a 
lenticular S0(6/7) galaxy, as indicated by its exponential brightness 
profile observed in several filters. Based on long-slit observations, 
these authors found different radial velocities for NGC 6908 lines 
in absorption (3113 ± 73 km s~^) and in emission (3060 ±16) 
km.s^^. Comparing with the radial velocity that they estimated for 
NGC 6907 (3190 ± 5 km s"^), the velocity difference between the 
two systems would be in the interval of 77 km s~^ to 130 km s~^. 
Summing the evidences of morphological asymmetry in the arm 
structure of NGC 6907, the behavior of its dust lanes and the tail 
of debris observed at low surface brightness, they conclude that the 
two galaxies are certainly interacting. 

Despite the convincing arguments in favor of the interacting 
nature of the two galaxies, one could still question whether the ab- 
sorption and emission lines observed in the direction of NGC 6908 
are indeed from that galaxy or originate from the rotating disk of 
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NGC 6907. Studies of rotation curves of spiral galaxies based on 
ionized gas (emission lines) and on the stellar component (absorp- 
tion lines) have shown that, for most galaxies, the rotation of the 
two components is similar, but local differences of velocities up 
to 100 km s"^ can occur (Saito et al. 1984; Fillmore et al. 1986; 
Kormendy & Westpfahl 1989; Vega Beltran et al. 2001; Pizzella 
et al. 2004 and Rhee et al. 2004). On the other hand, only lumi- 
nous SO galaxies are expected to present detectable emission lines. 
Since NGC 6908 is a low luminosity S0(6/7) galaxy (Ms=-17.4 , 
Madore et al. 2007) superimposed on the spiral arm of NGC 6907, 
the lines used to measure its velocity could well originate in the 
rotating material of the disk of NGC 6907. 

We present in this paper new observations, at radio and opti- 
cal wavelengths, that allow us to definitely settle this question, and 
to propose a scenario that describes in more detail the recent colli- 
sion between the two galaxies. We performed high resolution radio 
observations at 2 1 cm with the Giant Meterwave Radio Telescope 
(GMRT), and a detailed spectroscopic study over the field of the 
system NGC 6907 and 6908, with the Gemini North telescope. 

At 21 cm, using the aperture synthesis technique, the H I in- 
tensity maps and the velocity fields of both objects were obtained. 
The emission of the two objects could be distinguished and their 
velocities precisely measured. Furthermore the excess in velocity 
in the velocity field was used to estimate the time passed from the 
start of the interaction. 

The optical multi-object spectroscopy, employing the GMOS 
instrument, allowed the study of the emission and absorption lines 
of selected H II regions in the NGC 6907 disk and the nucleus of 
both galaxies. Each observed spectrum could be compared with its 
radio counterpart and the differences in absorption and emission 
velocities in direction of NGC 6908 could be explained by the idea 
that the emission lines were generated on the NGC 6907 disk and 
the absorption lines were formed on NGC 6908. 



2 THE FIELD OF NGC 6907/8 

Inside a field of 24' x 24' centered on NGC 6907 (Figure [B, only 
two extragalactic objects dominate the optical emission: NGC 6907 
itself, an SBbc spiral galaxy (Eskridge et al., 2002) with an absolute 
B magnitude Ms = -21.34 and a surface brightness of 22.29 mag 
arcsec"^ at half B total light (Lauberts & Valentijn 1989) and NGC 
6908, a S0(6/7) galaxy with Mb=-17.4 (Madore et al. 2007), 43.6 
arcsec north-west from NGC 6907. Garcia (1993) mentioned NGC 
6907 as the main member of a sparse group of 3 galaxies (IC 4999 
and IC 5005, without considering NGC 6908), also identified in the 
catalog of galaxy groups of Maia et al (1989), but these objects can- 
not be seen in the field shown in the figure. In the radio continuum 
at 1.4 GHz, NGC 6907 is also the brightest object, as can be seen 
from the NRAO VLA Sky Survey (Condon et al. 1998). The main 
properties of NGC 6907/8 galaxies are listed in Table[T] 



3 GEMINI OBSERVATIONS 

Optical observations were conducted at the Gemini-North tele- 
scope using the imaging and multi-slit spectroscopic capabilities 
of the Gemini Multi-Object Spectrograph (GMOS) in queue mode. 

Before the observations, images of NGC 6907 publicly avail- 
able were used to pre-select the best H II regions candidates on it as 
well to determine the optimum position angle for the observation 
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Figure 1. Optical image of the field centered in NGC 6907 galaxy and cov- 
ering 24 X 24 arcsec^ . It is overlaid by the GMRT continuum observation. 
The contours are associated with the levels 2, 4, 11, 18, 25, 33, 40, 47, 54 
and 61 mjy/beam. The external thicker contour con'esponds to the limit of 
detection of 2 sigmas of the noise level for the combined continuum im- 
age (~ 0.6 mJy/beam). Optical isophotal contour observed at the Gemini 
Telescope (gray) and elliptical fit (dashed line) are presented too. 

(243° or 63° - from north to east). This angle minimize the num- 
ber of selected objects that would have overlapping spectra on the 
GMOS detector. 

As frequently happens, the MOS observations were broken in 
two stages: one for imaging and the other for spectroscopy. The im- 
ages obtained in the first stage were employed to design the metal- 
lic mask with slits to perform the spectroscopy of the field and then 
used for photometric and morphological studies. 

As required in the queue mode a set of minimum conditions 
under which the observation could be executed were solicited. A 
minimum image quality of 70% was requested and for this rea- 
son observations at 2x2 binned mode were chosen, shortening 
the readout time. As a compromise between spectral resolution 
and amount of light detected all the slits were fixed with 1 arcsec 
widths. A standard slit height of 5 arcsec was chosen to improve 
the spatial sampling of the sky, since the typical source size of the 
candidates is almost 2 arcsec. 

All the observations were executed in better conditions than 
required, ensuring full photometric conditions in all images. The 
seeing was 0.7 arcsec and the sky brightness was as good as 21.7 
mag arcsec"^ in the g' filter and 20.1 mag arcsec^^ for the r' filter 

For all the observations fast tip-tilt guiding was per- 
formed using wavefront sensors for image motion compensation. 
Only one on-instrument wavefront sensor (OIWFS) guide star 
(GSC0691 101076) was needed for both imaging and spectroscopy 
observations. 

Optical imaging for the NGC 6907 field was executed on Au- 
gust 1st 2005 using the filters g' (g_G0301) and r' (r_G0303). To 
guarantee a continuous coverage of the NGC 6907 field we applied 
different dithering patterns to integrate a single frame in each filter 

The image files prepared in the first stage of the observations 



Table 1. Main properties of tlie galaxies NGC 6907 and 6908. The index 
opt is for the optical observations (combining the g' and r' filters) and rad 
for the radio observations. RA and Dec are the equatorial coordinates, 
Vb is the systemic velocity, a is the size, i the inclination, PA the posi- 
tion angle, fjL b the surface brightness. Type the morphological type, D the 
galactocentric Hubble flow distance, Shi the 21-cm Flux density, Ab the 
Galactic extinction in B-band, Vmax the maximum velocity rotation, r'-g' 
the color in these filters, B the absolute B-band magnitude. Lb the total 
B luminosity, M^y„ the dynamical mass and Mhi/I^b is the H i mass to 
B lununosity ratio. The number for the references are: (1) Eskridge et al. 
(2002), (2) Madore et al. (2007) and (3) Amores & Lepine (2005). 



Parameters 


NGC 6907 


NGC 6908 


RA (J2000) 


20''25"'06^63 


20''25™08^.97 


Dec (J2000) 


-24°48'33".0 


-24°48'04".3 


Vb [km/s] 


3182.4±3.9 


3147±21 


aopt [arcsec] 


90.0±1.2 


10.3±0.6 


a-rad [arcsec] 


211.8±1.9 




iopt [°] 


31.1±3.0 


74.7±3.1 


'^rad [ ] 


46.4±3.1 




PAopt [°] 


76.4±4.1 


2.8±1.4 


PArad n 


52.3±11.5 




Ms 


22.36±0.05 


20.81±0.08 


Type 


SBbc (1) 


S0(6/7)(2) 


D [Mpc] 


44.5±3.1 


43±16 


Shi [mJy.km/s] 


48.3±1.8 




Ab 


0.12±0.07 (3) 


0.12±0.07 


Vmax [km/s] 


215.8±6.7 




r'-g' 


-1.55±0.04 


-1.50±0.04 


B 


-21.69±0.13 


-17.43±0.15 


Lb [Lq] 


(6.9±1.3)xl0i" 


(6. 1 ±0.9) X 10** 


Mhi [Mq] 


(8.3±0.4)xl03 




Mdyn [Mq] 


(3.3±0.4)xl0" 




Mhi/Lb [Mq/Lq] 


0.12±0.01 





were used to select 33 H II regions candidates over the NGC 6907 
field, according to their brightness and geometry (Figure^. Stan- 
dard procedures using the GMOS Mask Making Software were em- 
ployed to design a mask with slits over each one of these candi- 
dates. A special slit was reserved for NGC 6908 (Slit 13). 

The spectroscopic observations were carried out in Septem- 
ber 2nd 2005 in a single exposure of 960 seconds using the 
B600_G5303 grating. This disperser and the slit width of 1 arc- 
sec guarantees a spectral resolution of ~5 A. No spectral dithering 
was used to cover the gaps on the final spectra, so special care was 
taken in choosing the central wavelength of 5000 A by consider- 
ing the systemic velocity of NGC 6907 and its expected rotation 
curve to avoid signal losses. The final spectra covered wavelengths 
between 3500 A - 6300 A. A CuAr arc lamp and a flat-field frame 
were observed in the sequence of the observation, while the tele- 
scope was still following the source, to guarantee the quality of the 
calibration frames, by avoiding any effect caused by torques over 
the instrument. 



3.1 Data reduction procedures for Gemini Observations 

There are great differences for imaging and spectroscopy GMOS 
reduction, but for both cases the Gemini Data Reduction Software 
was used. 

The reduction for images in filter g' and r' followed the same 
procedures. Individual dithered images of the field were corrected 
by overscan level, differences in gain for the detectors, bias and 
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Figure 2. GMOS image of NGC 6907 in the g' filter overlapped by the 
representation of the slits. The slit 13 was used to sample NGC 6908 and 
the slit 19 to sample the optical center of NGC 6907. Special squared slits 
signed with (A) were used to guarantee the alignment of the slits over the 
scientific targets. 

flatfield. After these steps images from different detectors were mo- 
saiced and all exposures were co-added in a clean image. 

Spectroscopic data reduction requires some extra cares. The 
first step was to prepare combined and normalized flatfield images, 
correcting this frames by overscan and bias. Since there was offsets 
between flatfield and the acquisition frames the normalized flatfield 
was used for finding the slit edges using the gradient method. For 
this observation, an offset of 6 pixels in the spatial direction was 
applied. 

Considering that a single frame was acquired in the spectro- 
scopic observations, the cosmic rays could not be removed by com- 
bining frames. In this case the task GSCRREJ was used. Residual 
effects on the spectral lines were corrected reprocessing only the 
slits affected by such correction, making the parameters for the cos- 
mic ray removing more flexible. 

An interactive process for the identification of spectral lines 
was set in each one of the spectrum analyzed. Multiple solutions 
were obtained for the same spectrum using the spatial distribution 
of the spectral lines, improving the wavelength calibration. 

For each slit a local sky spectrum was extracted to be used 
to subtract the sky lines and continuum from the target spectrum. 
Standard sky lines were used to verify if there was raw systemic 
error in the observation. 

After the visual inspection, the final spectra were generated by 
considering the on source integrated emission in the spatial direc- 
tion. 



4 GMRT OBSERVATIONS 

The interferometric mode of the Giant Meterwave Radio Telescope 
(GMRT) was used for 21 cm imaging and spectroscopy. A more 
detailed description about this instrument is given by Swarup et al. 
(1991) and Ananthakrishnan & Pramesh Rao (2002). 
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Table 2. Observational parameters for the GMRT observation. 



Parameter 



Value 



Date of Observation 
Pointing center (RA J2000) 
Pointing center (DEC J2000) 
Observing duration[hrs] 
Total on source time [hrs] 
Central frequency [MHz] 
Bandwidth per IF [MHz] 
Number of spectral channels 
Polarizations 

Frequency resolution [kHz] 
Velocity resolution [km.s^^] 
Flux calibrators 
Phase calibrator 
Bandpass calibrator 



2006 July 13 
20'"25™06^652 
-24°48'33".492 
4.7 
2.6 
1405.40 
4.0 
128 
2 

31.3 
6.7 

3C48 and 3C286 
1923-210 
3C48, 3C286 and 1923-210 



The observations were carried out with a total baseband band- 
width of 4.0 IVIHz divided into 128 channels centered at 1405.40 
MHz to cover the H I heliocentric velocity range of 2755 - 3610 
km s~^ (optical definition). A summary of the main observational 
parameters are presented in Table |2] 

The spectral resolution of the observation was 6.7 km s^^ per 
chatmel. Topocentric to heliocentric conversion of frequency was 
done with the help of DOPSET task to take care of the earth's mo- 
tion at the date and time of observation. 

The radio sources 3C286 and 3C48 were observed at the be- 
ginning and the end of the observation run as flux calibrator for 
~ 1 1 minutes each. VLA calibrator source 1923—210, which is an 
unresolved point source for GIVIRT (flux density of 2.0 Jy at 20 cm) 
was used as the phase calibrator. It was observed alternately for ~ 
5 minutes every ~ 20 minutes observing NGC 6907. 



4.1 Data reduction procedures for GIVIRT Observations 

The radio data were reduced following standard calibration and 
imaging methods employing the classical AIPS. Flux densities at 
the observing frequency were estimated using the standard primary 
calibrators 3C286 and 3C48, also used to evaluate the flux densi- 
ties and phases of the secondary calibrator. For 3C286 the estimated 
flux density was 14.81 Jy and for 3C48 it was 16.03 Jy. Taking into 
account the results of Omar et al.(2002) for 3C48, the accuracy in 
the flux calibration must be better than 3%. Since there was not any 
spectral features for 3C48, 3C286 and 1923—210 at the observing 
frequency, all of them were used to bandpass calibration. 

Standard spectral line data reduction procedure was followed 
to remove bad data, to calibrate the gains and bandpass to produce 
the final calibrated visibilities for the target field. 

A preliminary datacube was generated at spatial frequencies 
below 5 kA to verify which channels would contain H I emission. 
Visual inspection, channel by channel, enabled to register the de- 
tection of extended H I emission within the channels 41-105. 

The polyhedron imaging technique was used to generate con- 
tinuum images by averaging 42 line-free channels (ranging from 7 
to 36 and 112 to 125), with data self-calibrated in both phase and 
amplitude. Spatial frequencies below 5fcA were used and a Gaus- 
sian taper function was applied to weight down long baselines, re- 
ducing the sidelobes. In this case the overall noise level was 2 mjy. 
The continuum flux density was subtracted from each individual 
channels and the final self-calibrated solutions were applied to the 



Table 3. Parameters derived from the elliptical fitting to optical and radio 
isophotal data at 3ct of the sky or noise level; a and b are the major and 
minor axis of the ellipse, in arcsec, qq and Sq are the equatorial coordinates 
for the center (J2000), i is the inclination, PA is the position angle, both in 
degrees. 



Para- 


NGC 6907 


NGC 6907 


NGC 6908 


meter 


Optical (r'+g') 


Radio (5kA) 


Optical {r'+g') 


a 


88.2±1.8 


212±19 


10.5±0.8 


b 


77.4±1.2 


148±10 


3.2±0.2 




20'^25'"06=.29 


20''25'"08^70 


20''25'"08''.96 


Sq 


-24°48'16".3 


-24°48'43".6 


-24°48'03".4 


i 


29.1±4.2 


46.4±3.1 


74.1±4.5 


PA 


74.1±6.6 


52.3±11.5 


3.1±2.1 



line-data to make an image cube with a resolution of 52" x 39". 
To enhance the sensitivity to extended features the data points were 
"natural-weighted" and the resulting spectral cubes were cleaned 
for signals greater than 2.5 times the rms noise in the channel im- 
ages. 

To improve the spatial resolution for the analysis, several it- 
erations of self- calibration were executed to make image cubes at 
different UV-ranges adapting the previous procedure. 

The relationship between frequency and spectral-line velocity 
was set by means of the standard task ALTSWTCH. IVIoment maps 
were made from the resulting data cubes by using the AIPS task 
IVIOIVINT (see the contours in Figures |6] and |9] for the 0th and 1st 
moment map respectively). 



5 RESULTS 

5.1 Optical Photometry 

The photometric calibration for the images in the filters g' 
and r' was made using the photometric standard stars in the 
IVIARK-A field (Landolt 1992), supplied by the Gemini base- 
line (gn-cal2005061 1) and the usual procedures of the DAOPHOT 
IRAF tasks. Solutions including the atmospheric extinction were 
considered to determine the zero point magnitudes for the g' (mo = 
31.72) and r ' (mo = 31.60). The results are in excellent agreement 
with the median atmospheric extinction coefficients for IVIauna 
Kea, with differences lower than 0.4% 

The limit of detection of the galaxies was taken as the 3cr limit 
of dispersion of the local sky level. For NGC 6907 these limits are 
24.7 ± 0.2 mag arcsec"^ (g') and 23.6 ± 0.2 mag arcsec"^ (r'). 
For NGC 6908 they are 22.6 ± 0.3 mag arcsec"^ and 21.4 ± 0.3 
mag arcsec"^ respectively to the same filters. 

Elliptical fittings in these isophotal limits, masked from the 
contribution of stellar fields, the interarm regions and the disturbed 
parts of the galaxy provided the projection parameters of NGC 
6907 and 6908 (Table[3). 

The integration of all flux inside these isophotes in each fil- 
ter gives the g' and r' magnitudes. These values were converted 
to the standard B-magnitude using the expressions by Fukugita 
et al. (1996) and Hook et al. (2004), corrected by the foreground 
Galactic extinction (Amores & Lepine 2005) and masked from the 
foreground contributions. The results were a total B-magnitude of 
11.55 ± 0.03 and a mean effective surface brightness of 22.36 ± 
0.05 mag arcsec " ^ for NGC 6907 and a total B -magnitude of 1 5 . 8 1 
±0.15 and a mean effective surface brightness of 20.8 1 ± 0.08 mag 
arcsec"^ for NGC 6908. 
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Figure 3. g' and r' radial surface brightness profiles. Special limits men- 
tioned in the text are indicated by vertical or horizontal lines over the graph. 



The radial surface brightness profiles were obtained perform- 
ing the integration of the flux on rings with the same size of the 
seeing (0.7 arcsec) and following the parameters of projection ob- 
served for the galaxy in each filter. For NGC 6907 (Figure |3j the 
central variation of the brightness is compatible with a small spi- 
ral structure (usually identified as a bar) within (7.4 ±2.1) arcsec. 
The exponential behavior is only affected at radius near NGC 6908 
emission, returning back to the exponential behavior for larger dis- 
tances. For NGC 6908 the surface brightness profiles are compat- 
ible with Madore et al. 2007, but the extended central emission of 
NGC 6908 suggests the presence of a composed nucleus. 



5.2 Optical Spectroscopy 

For the spectrophotometric calibration the sensitivity function was 
achieved using the flux standard star G191B2B supplied by the 
Gemini baseline, according to the specifications for the NGC 6907 
spectral range coverage. 

The IRAF package RVSAO was used to identify the spectral 
lines observed in the H II regions of NGC 6907 by cross-correlating 
a template spectrum of this kind of object to the observed spectra. A 
set of shifts and errors were determined for each spectral line and 
the radial velocities of all lines in each spectrum were combined 
by discarding the weakest lines and weighting the velocities with 
their errors and equivalent widths. The results of these process are 
summarized in Table |4] 

Special care was taken with the NGC 6908 spectrum (Figure 
|4]l, since its optical depth is unknown and it can have contributions 
coming from NGC 6907 in its own spectrum. Taking into accoimt 
that NGC 6908 is a lenticular SO galaxy, then absorption lines must 
provide the most reliable information about its contents. Visual in- 
spection was done to recognize spectral absorption lines and in- 
dividual fitting profiles were performed considering Gaussian and 
Lorentzian models to obtain the central wavelength of the lines (see 
Table|5j. Emission lines were detected inside the absorption lines, 
specially in the hydrogen lines (Figure [5j, confirming the spectra 
superposition. 



Table 4. Observed radial velocities in direction of the selected H II regions 
in the optical (Vopt) and in 21cm (yni)- The ID is the slit number, RA and 
Dec are the equatorial coordinates (J2000), r is the de-projected distance 
from the center, considering the projection parameters of the radio for which 
the uncertainties are lower than 0.2 kpc (Ho = 73 ± 5) km s^^/Mpc. 
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Figure 4. NGC 6908 integrated spectrum, smoothed by removing Fourier 
components of higher frequencies. Absorption and emission lines are iden- 
tified by the arrows. 



5.3 Continuum 

A more detailed view of the continuum emission of NGC 6907 
shows a north-east elongated structure, with an asymmetric emis- 
sion which does not provide good elliptical fitting (Figure [TJ. The 
interaction between NGC 6907 and 6908 must be the main reason 
for such asymmetry. Almost all the continuum emission is confined 
to an isophote of 1.9 mjy/beam with approximately the same size 
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Figure 5. Detail of each absorption line (continuous line) and the fitted pro- 
file (dashed line). A Gaussian or Lorentzian profile was chosen depending 

2 

on which minimizes the x P^r degree of freedom. 



Table 5. Central wavelengths and velocities for the absolution specti'al lines 
observed in direction of NGC 6908. Aq is the rest wavelength. \abs the 
peak centroid, and v^hs is the velocity associated to displacement of the 
observed spectral line. 



Absorption 


Au 


^abs 




Line 


[A] 


[A] 


[km/s] 


K 


3933.70 


3975.91±0.36 


3217±13 


H 


3968.50 


4011.31±0.28 


3234±16 


H<5 


4101.70 


4145.49±0.24 


3201±28 


H7 


4340.50 


4386.68±0.31 


3189±28 


H/3 


4861.33 


4914.56±0.54 


3283±21 


Mg 


5175.36 


5224.40±0.28 


2841±17 


CaFe 


5268.98 


5323.46±0.49 


3099±21 


Na 


5892.50 


5954.66±0.25 


3162±33 



than the optical emission centered at ac=20''25"'07^.28 and 5^=- 
24°48'33".9 . The total continuum flux measured for NGC 6907 
is 112 mjy, in absolute agreement with Condon et al. (1998) and 
Vollmer et al. (2005), if the same isophotal limit is considered. The 
radial profile of this emission exhibits the expected behavior for an 
exponential disk inside a de-projected radius of 72 arcsec. 




lOs.O 20h25iii05s.O 
Right Ascension (J2000) 



Figure 6. H I column density map over the g' image of NGC 6907 and 
6908. The contour levels are at 0.3, 8.8, 17.3, 25.9, 34.4, 42.9, 51.4, 60.0, 
68 .5 and 77.0 in units of lO^^ cm" 2. 



5.4 H I Emission 

Subtracting the continuum image from the image cubes composed 
for UV-ranges up to 40 kA it was possible to recover the spatial 
H I distribution for each channel, resolving structures larger than 
lOx 12 arcsec. That is enough to sample the integrated contribution 
of the components inside the galaxy NGC 6907. 

The bidimensional integration of the H I emission detected 
above 2a of the noise level in all channels composes the 0th mo- 
ment. In Figure |6] this map was converted to the H I column den- 
sity assuming the relation by Spitzer (1978). It reveals a substan- 
tially more extended spiral structure than observed in optical im- 
ages and the emission has two extended peaks coincident with the 
spiral arms. 

Elliptical fitting for the isophote associated with a column den- 
sity of 2.8 X 10^^ cm~^ provides the extension of the radio emission 
as (3.5 ± 0.3) arcmin or almost 2.4 times larger than the optical 
emission (Table[3}- Neither the optical center nor the radio center 
are coincident. 

The integration of the column density on elliptical rings along 
NGC 6907 provides the radial column density profile presented in 
Figure [7] As expected, there is depletion of H I on the central re- 
gions of NGC 6907, since the molecular hydrogen is formed in a 
denser environments. After the peak emission at 10 kpc an expo- 
nential fall occurs up to 28 kpc, where the column density rises up 
to 34 kpc and falls back again. Such distribution is quite common 
as can be seen in Broeils & van Woerden (1994), and the minimum 
inside the galactic disk can be related to the gas redistribution in 
the corotation radius. 

On the other hand, the integration of all the emission for each 
channel of the datacube after continuum subtraction provides the 
global H I emission profile for the galaxies NGC 6907 and 6908 as 
it can be seen in Figure [8] 

The double peaked emission is expected for rotating disks. 
Nevertheless, the different levels of the peaks indicates the pres- 
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Figure 7. Radial H I column density profile. The open circle with error bars 
informs the mean enor in data. 
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Figure 8. Global H I emission profile of NGC 6907 and 6908. The bottom 
dotted lines limit the noise level while the top dotted line indicates the max- 
imum emission h. Inside the distribution two horizontal dashed lines are 
used to mark the H I Hne- widths estimated at 20% (w2o) and at 50% (W50) 
of the peak intensities; v is the central weighted mean velocity distribution. 



ence of other components in the distribution of velocities measured 
in the field. These effects can be attributed to the spatial velocity of 
NGC 6908 projected on the line-of-sight and to the intrinsic asym- 
metries on the velocity field of NGC 6907. 

Richter & Huchtmeier (1991), Chengalur et al (1993) and 
Doyle et al (2005) observed these galaxies in lower resolution with 
different single dishes, what can possibly account for the main dif- 
ferences in the registered line profiles. The mean differences are 
lower than 3%, including the flux density (48.3 ±1.8 mjy.km s~^), 
suggesting that a beam smearing is the main responsible for the dif- 
ferences found at the flux peak. 

The weighted mean velocity associated with the H I profile 
(3182.4 ± 3.9 km s^^) is absolutely compatible with the velocity 
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Figure 9. NGC 6907 velocity field superimposed over its image in g'. 
Dashed lines represent the part of the velocity field going away from the ob- 
server while the opposite is represented by the continuous lines. The thicker 
crossing lines represent the kinematical axis. The velocity field was gener- 
ated using an UV range of 5fcA and a resolution of 52" X 39", to balance 
the sensitivity to the diffuse emission and the spatial resolution. 



measured by the robust mean of the entire 1st moment map (3183.2 
± 6.7 km s~^). Taking into account the Hubble flow distance (sup- 
posing H = (73 ± 5) km/sec/Mpc), the galactocentric distance of 
this galaxy is estimated as (44.5 ± 3.1) Mpc. Using this distance a 
H I mass of (8.3 ± 0.4) x lO** M© was calculated. 

Considering the distribution of velocities associated with the 
bidimensional emission in each frequency of the 1st moment map, 
then the observed velocity field of the galaxy can be extracted (Fig- 
ure |9ll. 

One can note that the velocity field is strongly oriented by 
the spiral arm in the internal regions. In fact, the high density of 
the H I emission in the spiral arms (from 10 to 100 times higher 
than the interarm regions) could generate a bias in the observed 
velocity field, if the velocity from the disk does not correspond to 
that measured in the spiral arms. 

The kinematic major axis is not perpendicular to the kinematic 
minor axis (Figure [9}. That can be easily verified by comparing 
the linear fit over the coordinates of the points with the systemic 
velocity (minor axis) and the coordinates of the nearest points to 
the minor axis of all isovelocity levels, which samples the major 
axis in an axisymetric rotation disk. The angle between these axis 
is (78.2 ± 3.1)° and the point where they cross each other can be 
used to define the kinematic center (0:^=20'' 25" 06''. 40 and Sk=- 
24°48'38".3). 

Since the the global H I profile suggest that a rotating disk 
must be the main component of the movement of the system, then 
the rotation curve of the galaxy NGC 6907 can be extracted sup- 
posing that the elliptical fitting made on the H I gas distribution 
represent the projection of the NGC 6907 disk on the sky and its 
gas follows circular orbits. In this situation the observed line-of- 
sight velocity ignoring vertical and radial movements would be: 



^(q, S) = Vo+ Vrot{r, 9) sin(j) cos(6l) 



(1) 
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Figure 10. Rotation curves extracted for the receding and approaching re- 
gions of the NGC 6907 velocity field. The hatched area correspond to the 
positions occupied by NGC 6908 in the receding part of NGC 6907 and two 
vertical dashed lines represent elliptical fit limits in g' and r' . For the stellar 
rotation curve it was supposed a mass-to-light ratio of 1.5. 



taking into account that: 



^Qg^gi) ^ -(q - Qo) sin(0) cos(^) + {5- 5q) cos((/.) 



sm{9) 



-[a — ao) cos((/)) cos(5) — (5 — So) sm{(j>) 



(3) 



where (ao, 5o) is the kinematic center, Vo is the systemic ve- 
locity, 4> is the position angle of the receding semi-major axis, i is 
the inchnation between the normal to the plane of the galaxy and 
the line-of-sight and r and 6 are the polar coordinates inside the 
plane of the galactic disk. 

Note that we introduce the term cos((S) multiplying all dis- 
tances measured in right ascension (what it is not usually made in 
the literature), otherwise it would overestimate the galactocentric 
radius, depending on the declination of the source and the position 
angle of each target relative to the center of the galaxy. For the NGC 
6907 these errors are predominantly of the order of 12%. 

Using the combination of the velocity fields made with maps 
covering UV-ranges from IkA to 40kA, taking into account that 
the most reliable information about the rotation curve is extracted 
near the kinematic major axis ±30°), then the previous equa- 
tions can be used to determine Vrot, r and 6 for the receding and 
approaching regions of the galaxy (see Figure [TOt. The tilted ring 
model was applied to the NGC 6907 velocity field, considering the 
same limits used in the previous method. Variations of the incli- 
nations and the position angles with the distance are lower than 
the uncertainties for the position angle and inclination of the whole 
galaxy so its results can be ignored. A raw estimate for the dynam- 
ical mass of NGC 6907, considering the last point in the rotation 
curve is (3.3 ±0.9) X 10" Mq. 

Comparing the radial velocities from optical and H I spectra 
observed in the same directions (Figure [TT] and Table |4] it is pos- 
sible to verify the overall correspondence between the optical and 
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Figure 11. Comparison between optical radial velocities observed with 
GMOS and radio radial velocities observed with the GMRT. The straight 
continuous line represent the linear least squares regression for the observed 
regions with azimuthal angles near the major axis (AO < 30°). 



radio rotation curves {vopt = (1.00 ± 0.08)vhi ± (1 ± 240)) and 
verify that the radio projection parameters provides the best corre- 
spondences. The high dispersion of individual points is related to 
the higher uncertainties to extract radial velocities azimuthaly far 
away from the major axis. In spite of that, the mean dispersion is 
inferior to the limit of 50 km s~^, refereed by Bosma (1981) as 
the limit related to large-scale asymmetries of velocity fields. The 
exclusion of the objects in larger azimuthal angles relative to the 
major axis {A9 = 30°) provides almost the same fit. In particular, 
the point representing the slit which sampled NGC 6908 presents a 
higher velocity dispersion than the others, possibly due to the spa- 
tial motion of NGC 6908. 



6 DISCUSSION 

A first view on the NGC 6907 velocity field (Figure|9) reveals that 
NGC 6908 is over the approaching part of it. The extraction of the 
rotation curve in this part of the galaxy reveals contributions of ra- 
dial components of the NGC 6908 movement, what does not occurs 
in the receding part of the galaxy. This verification explains the dif- 
ferent peak levels in the global H I emission profile shown in Figure 
[8] and in the rotation curves presented in Figure [TO] Taking into ac- 
count that the H I global profile preserves the behavior expected for 
a rotating disk and different emission profiles agree with a exponen- 
tial disk, then the effects of the interaction on the velocities can be 
considered as local. Consequently the rotation curve extracted from 
the receding part of the velocity field is the real (undisturbed) NGC 
6907 rotation curve. 

Comparing the raw observed rotation curves from the receding 
and approaching parts of NGC 6907 in direction of NGC 6908, it 
is possible to see that the relative emission velocity between these 
galaxies (130 ± 17 km s^^, as measured by Madore et al. 2007) 
, coincides exactly with the expected velocity for the undisturbed 
rotation curve, and not with the velocity measured in radio for the 
approaching velocity field (where NGC 6908 can be found). This 
apparent contradiction can be easily understood considering that 
the emission line velocity measured for NGC 6908 is in fact from 
the rotating material in the NGC 6907 disk. NGC 6908, as a lentic- 



The pair of galaxies NGC 6907 and 6908 9 



ular galaxy, has less gas and new hot stars to heat it. So the emission 
lines used by Madore to measure the NGC 6908 velocity only could 
come from the material in the NGC 6907 disk. 

The idea that the emission lines come from the NGC 6907 is 
corroborated by the spectrum extracted in direction of NGC 6908 
(Figure |4]l. It has a strong continuum superimposed by few emis- 
sion lines, for which the radial velocity is (3092 ± 52) km s~^, and 
by several absorption lines, associated to a radial velocity of (3147 
± 21) km s^^ (Table[5). Relative velocities spectroscopically mea- 
sured in relation to the NGC 6907 center (3209 ± 10 km s~^) are 
117 ± 53 km s"^ for emission lines (compatible with the NGC 
6907 rotation curve in that direction) and (62 ± 23) km s~^ for 
absorption lines, in perfect agreement with Madore et al. (2007). 

Another confirmation about the origin of the emission lines 
is that the broaden hydrogen absorption lines (typical of lenticular 
galaxies) contain vestiges of narrow emission lines inside it (Fig- 
ure |5](. These emission lines have the same radial velocity of the 
rotation curve of NGC 6907 at that radius. That is the reason why 
different references and observations do not agree with the presence 
of certain hydrogen lines (they can be hidden in lower resolutions 
observations) and why Hf3 is observed in absorption while Ha in 
emission. Since Ha line is more intense than Hf3 the equivalent 
absorption line formed in a different source is not deep enough to 
hide it. 

Verifying the distribution of H I emission as a function of the 
velocity in the direction of NGC 6908, it is possible to distinguish 
three tied Gaussian components (Figure[T2t. Considering the rela- 
tive velocities measured with respect the radio systemic velocity 
expected for the nucleus (since the spectrum extracted in direc- 
tion of the slit 19 does not sample the radio nucleus), the centroid 
of each component can be interpreted as associated with different 
structures in the interaction between NGC 6907 and 6908. The first 
one, at 3120 km s^^, is associated with NGC 6908, since it has 
the same velocity of the absorption lines generated by its stellar 
content, the second one, at 3065 km s^^, is produced by the ex- 
cited material from the disk of NGC 6907 during the interaction 
with NGC 6908 and the third one (main component) related to the 
higher velocity gas left behind NGC 6908 due the interaction with 
NGC 6907. 



6.1 A scenario for the collision between the two galaxies 

Using the velocity contours above the approaching velocity of 220 
km s~^, it is possible to trace the disturbed regions of the velocity 
field and see a converging distribution of isophotes over NGC 6908, 
whose gradient is higher in direction of the rotation of NGC 6907, 
following the spiral arms (Figure fTSl. 

The striking triangular shape of the isophotes associated with 
high-velocity gas could be explained by the following scenario. 
NGC 6908 crossed the disk of NGC 6907 in a recent past, and very 
probably the tidal forces created a bridge of gas between the two 
galaxies, similar to the Magellanic stream in our galaxy. The gas 
left behind NGC 6908 disk by the collision explains the excess of 
gas approaching us. In our galaxy, at the solar Galactic radius, the 
time required for a gas cloud launched in the direction perpendicu- 
lar to the plane to revert its motion and fall back onto it is about 35 
Myr (see eg. Lepine & Duvert, 1994). This can be taken as an order 
of magnitude for the present case. Let us now consider the velocity 
component of NGC 6908 in the plane of the sky. This component is 
directed approximately towards the center of NGC 6907, as we can 
infer from the fact that presently NGC 6908 is situated closer to the 
center of NGC 6907 than the marks left in the gas of the disk by 
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Figure 12. H I spectrum in direction of NGC 6908 (gray stepped curve). 
The dashed vertical lines mark special velocities indicated in the picture. 
The horizontal dotted lines limit the noise level. Three components were 
fitted for this spectrum: the first associated with NGC 6908 by a dashed line, 
the second related to the NGC 6907 disk emission lines (dot-dashed line) 
and the third related to the gas left behind NGC 6908 due the interaction 
(dot-dot-dashed line). The sum of these components is represented by the 
smoothed black curve overlaying the observed H I profile (step gray line). 
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Figure 13. Contour of isovelocities higher than 220 km s~^ in the ap- 
proaching region of the NGC 6907 velocity field overlaying the g' image of 
this galaxy. The interval between contours is 5 km.s~^. 



the crossing event. The gas stream produced by the collision must 
also have a component of velocity in the same (almost radial) di- 
rection, which explains why the triangular isophotes point towards 
NGC 6908. Finally, let us analyze the component of velocity of the 
stream in the plane of the sky, but in the direction of rotation of 
NGC 6907. The stream can, in a simplified way, be considered as a 
mixture of gas with velocities ranging from that of NGC 6908 (with 
almost zero component in the direction of rotation of NGC 6907) 
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to the gas which resulted to have about the same rotation velocity 
of the disk of NGC 6907. This explains the elongated shape of the 
isophotes in the direction of rotation, and at the same time, allows 
us to estimate in a new way the time elapsed since the collision. The 
mean de-projected distance of the two extremes of the gas with per- 
turbed velocity in Figure [T3]is about (104.8 ± 17.9) arcsec, where 
the undisturbed rotation curve has the velocity of (200.0 ± 6.7) 
km . Assuming Hq = (73 ± 5) km s~^/Mpc, this distance 
is (22.6 ± 4.1) kpc. Since the angle between the two extremes is 
(17.7 ± 0.3)°, then the initial disturbed material would have trav- 
eled (7.0 ± 1.2) kpc azimuthally. From this we derive that the colli- 
sion occurred about (3.4±0.6) x 10^ years ago, which is consistent 
with order of magnitude previously estimated. 



7 CONCLUSIONS 

The combination of the GMRT radio observation and 
GMOS/Gemini optical observation provides a more complete 
view about the interaction between the galaxies NGC 6907 and 
6908. 

Photometric images observed at the Gemini Telescope and the 
H I maps observed at the GMRT reveals extended structures that 
can be understood as tidal debris in optical images, but resemble 
the continuation of the spiral arms in radio. By fitting ellipses to 
the last isophotes registered to each one of these observations for 
NGC 6907 we verified that the 21 -cm emission is 2.4 times more 
extended than the optical emission. The orientation of the fitted el- 
lipses provide different parameters of projection for optical and ra- 
dio observations, the differences being 15.3° for the inclination and 
24. 1° for the position angle. Geometric, photometric and kinematic 
centers are not coincident, reinforcing the evidences about the ef- 
fect of an interaction in NGC 6907. 

The spectra observed using the GMOS spectrograph are co- 
herent with those observed in 21 -cm in the same direction of the 
slits. The parameters of projection in radio provides the best cor- 
respondences between the optical and radio velocities. The details 
of the emission and absorption lines observed in direction of NGC 
6908 show that they have different origins. First, because the ob- 
served emission lines are not typical of lenticular galaxies like [Oil] 
(3727 A), [OIII] (4959 A and 5007 A) and the Hel (5876 A), spe- 
cially for those with low luminosity like NGC 6908. Second, be- 
cause the radial velocity observed in absorption is different from 
the velocity emission lines, which are coherent with the undisturbed 
rotation curve of NGC 6907 expected for that direction. Finally, 
because inside the hydrogen absorption lines it is possible to see 
traces of narrower hydrogen recombination lines with a Doppler 
effect suitable for the NGC 6907 rotating disk, confirming their 
different origin. 

The H I column density distribution presents a discontinuity in 
its radial profile which can be can be related to the gas redistribu- 
tion at the corotation radius. The details of the global H I emission 
profile reveals that the rotation is the main component of the H I 
velocity distribution, and its asymmetry indicates the contribution 
of the NGC 6908 motion is local and in velocities near to those 
found in the NGC 6907 rotating disk. The total H I mass obtained 
was (8.3 ±0.4) x 10^ A/©. 

Due to the dominant contribution to the H I emission in the 
spiral arms they can be the main contributors for the velocity field 
asymmetries, verified by the non perpendicularity of the kinematic 
axis (difference of 12°). 

The differences in the rotation curve extracted from the reced- 



ing and approaching parts of the velocity field can be interpreted as 
consequence of line-of-sight components produced by the interac- 
tion with NGC 6908. Assuming that the interaction only affected 
the side where NGC 6908 is observed in projection, and that the 
velocity field of the opposite side is the real rotation curve of NGC 
6907, we were able to get a raw estimative for the dynamical mass 
of (3.3 ±0.4) X 10" Mq. 

Looking at the H I spectrum in direction of NGC 6908 it is 
possible to distinguish three emission components: one coincident 
with the relative velocity for NGC 6908, another coincident with 
the excited gas in NGC 6907 disk and the last one in agreement 
with the higher relative velocities gas left behind NGC 6908 by dy- 
namic friction. Using the difference in the rotation curves extracted 
from the disturbed region of the velocity field and that one away 
from this region, we found the limits for the disturbed material on 
the NGC 6907. Supposing that this material has drifted inside the 
NGC 6907 disk, according the rotation curve, we estimated the time 
when the interaction started in (3.4 ± 0.6) x 10^ years. 

Another work is in preparation to discuss the metallicity distri- 
bution in the NGC 6907 and its possible cormection with the coro- 
tation radius. 
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